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We solve for the behaviour of cosmological perturbations in the Dvali-Gabadadze-Porrati 
(DGP) braneworld model using a new numerical method. Unlike some other approaches in 
the literature, our method uses no approximations other than linear theory and is valid on 
large scales. We examine the behaviour of late-universe density perturbations for both the 
' self-accelerating and normal branches of DGP cosmology. 
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1 Introduction 

The Dvali-Gabadadze-Porrati (DGP) ^2111 j^iodel postulates that we live in a 4-dimensional hy- 
' persurface in a 5-dimensional Minkowski bulk. General Relativity (GR) is recovered at small 

scales (smaller than the crossover scale r^) due to the inclusion of an induced gravity term in 
the action. 

""^^I This model has two distinct classes of cosmological solutions.^ One of them exhibits acceler- 

^ ■ ated expansion at late times without the need to include any exotic cosmological fluids, such as 

■ dark energy, or any brane tension that acts as an effective 4-dimensional cosmological constant. 
c/3 , Hence, this branch of solutions is called "self-accelerating" . To explain the observed acceleration 

we require Tc ~ H^^, where Ho is the current value of the Hubble parameter. It is expected 
that structure formation will help to distinguish the self-accelerating DGP universe from dark 
energy models based on 4-dimensional GR. This is because the growth of cosmological pertur- 
\jQ I bations is very sensitive to the existence of an extra dimension. A full 5-dimensional treatment 

0^ ■ is required to model these perturbations, which is why obtaining observational predictions for 

the behaviour of fluctuations in the DGP model is technically challenging. 

Several authors have considered the problem of the dynamics of perturbations in the DGP 
I model, but they have all relied on some sort of approximation or simplifying ansatz. One example 

• of this is the quasi-static approximation scheme, which solves the perturbative equations 

of motion by focussing on the extreme subhorizon regime. 

In this paper we present a complete numerical analysis of the evolution of scalar pertur- 
^ ' bations in the DGP model. Mathematically, the problem involves the solution of a partial 

■ differential equation in the bulk coupled to an ordinary differential equation on the brane. A 
numerical method for dealing with such systems has previously been developed for cosmologi- 
cal perturbations in the Randall-Sundrum (RS) model. However, the DGP problem is more 
complicated than the RS case due to a non-local boundary condition on the bulk field. Hence, 
the algorithm used in this paper represents a significant generalization of the one used in the 
previous case. 

As alluded to above, there is another "normal" branch of solutions in the DGP model. We 
cannot explain the late time accelerated expansion of the Universe using the normal branch 
without including an effective cosmological constant induced by the brane tension a. However, 
by allowing a non-zero a, normal b ranc h solutions can mimic dark energy models with the 
equation of state w smaller than — I.ISISI Unlike 4-dimensional models that realize w < —1 hy 
the introduction of a phantom field, the normal branch of DGP cosmology is ghost-free. This 
unique feature is what motivates us to numerically study the perturbations of the normal branch 
of DGP cosmology. 
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Based on the work^done in collaboration with Kazuya Koyama, Sanjeev S. Seahra and Fabio P. Silva. 



2 Background solution 



In the DGP model the brane dynamics is governed by 
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H 

a 2rc 



(1) 



^ = -3{l + w)pH, (2) 

where a{t) is the scale factor of the brane universe, normalized to unity today, a dot denotes 
the derivative with respect to the proper time along the brane, t, and Vc = kI/2k\. We assume 
that the stress-energy tensor of the brane matter is of the perfect fluid form. Note that, as in 
GR, the stress-energy tensor is conserved. The e = ±1 parameter reflects the fact that when we 
impose the Z2 symmetry across the brane, we have two choices for the half of the bulk manifold 
we discard. When e = +1, we have that Hr^ ~ 1 when the density of brane matter is small, 
|p + o"! ^ ^r~^. This implies a late-time accelerating universe, which is why the e = +1 case 
is called the self-accelerating branch and the e = — 1 case is called the normal branch. 

3 Master equations governing perturbations 

In this section we describe the formulae governing scalar perturbations in the late-time matter 
dominated universe. We take the matter content of the brane to be a dust fluid (i.e., cold dark 
matter), w = 0, therefore we have p oc . 

It can be shown that the dynamics of all the scalar perturbations of the bulk geometry can 
be derived from a single scalar bulk degree of freedom. After Fourier decomposition, we find 
that the mode amplitude 17 = Vl{u,v) of this master field obeys 

= ^_A^ + ^O, (3) 
dudv 2v du 4f 2 ' 

where u and v are dimensionless null coordinates. 

We also find the following boundary condition for Q: 

(9y^h - -^^b + ^f^b + A, (4) 

the following equation of motion for the density contrast of the cold dark matter. A: 

A + 2HA- ^4pj2^ = -'-^^h-, (5) 
and the following expressions for the two metric potentials ^ and ^: 
_ 4pa^-fi £71 ^ e(^^jf3FV)7i 

* = 2^^ + I^^^ - -^^"^ + 

where ilb = ^b(^) = ^(^b(0) ^^b(i)) and {dyO,)^^ are the values of the bulk field and its normal 
derivative, respectively, evaluated at the brane. In these expressions, the dimensionless 7-factors 
are functions of H and H. 

The bulk wave equation ([3]), boundary condition ([1]) and ([5]) are the equations we must solve. 
Once we know A and the metric perturbations $ and ^ can be obtained by differentiation. 
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Figure 1: Linear growth factor and alternate gravitational potentials ^± from simulations and the QS approxi- 
mation in the self-accelerating branch (left) and normal branch with Q.r^ — 0.05 (right). We normalize g{a) and 
"1>_ to unity at early times. For comparison, we also show the relevant results for the concordance ACDM model 

with = 0.26 and f^A = 0.74. 



4 Scalar perturbations in the self-accelerating universe 

In this section, we concentrate on the w = a = and e = +1 DGP scenario as a model for the 
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late-time accelerating universe. By examining probes of the expansion history, it has been 
found that 

^r. = j^ = 0.15±0m, (8) 

at 95% confidence. We use the best fit value for ^Ir^ in our simulations. 

For all plots in this paper, we select the bulk field to be zero and the brane field non-zero 
initially. We have also simulated several different choices of initial data, such as the bulk field 
being constant along the initial null hypersurface, and have found that the simulation results 
remain the same as long as the initial time is early enough. This is analogous to what happens 
in the RS case. ' 

In a previous work, ^ a 'quasi-static' (QS) approximation was developed to describe the 
behaviour of DGP perturbations whilst well inside the cosmological horizon, k ^ Ha, and with 
physical wavelengths much less than the crossover scale, a <C krc- Here, we compare the QS 
approximation to our simulations to determine just how large k must be for it to be valid. 

In Fig. [T] (left), we compare simulation results versus the QS approximation for the linear 
growth factor g{a) = A(a)/a and the alternate gravitational potentials $± = |(<I> ± ^f). We see 
that the simulation results are consistent with the QS approximation for k > /iMpc~^. On 
larger scales, the potential which determines the integrated Sach- Wolfe (ISW) effect, shows 
more suppression than the QS prediction. 



5 Scalar perturbations in the DGP normal branch 



We now turn our attention to the behaviour of density perturbations in the normal branch of the 
DGP model. Unlike the e = +1 case, this branch does not naturally have a late time accelerating 
phase. So, in order to be made consistent with observations, we must allow for the brane to have 
a nonzero tension that acts as an effective 4-dimensional cosmological constant (we call this the 
ADGP model). Assuming that the matter sector is CDM-dominated, the Friedmann equation 
for this scenario follows from the general form ([1]) with e = — 1 and w = 0. The background 

1 

dynamics has been compared to observations of H{z), and the following parameter values 
were found: ^ 

= ^ = 0.23 ± 0.04, = — ^ < 0.05, (9) 

at 95% confidence. Here, po is the present day CDM density. Note that the observationally 
preferred value of 0^^. is zero. Since the DGP model goes over to GR in this limit, this implies 
that ACDM gives a better fit to the data than ADGP. We will assume the best fit value of 0.23 
for rim. 

In Fig. [T] (right), we compare the results of our simulations to the QS approximation and 
ACDM in the case = 0.05. As in ^ we find that the simulation results are fairly insensitive 
to initial conditions provided that the initial data surface is set far enough into the past. In 
contrast to the self-accelerating case, we find that the linear growth factor and potential 
are generally larger than in the ACDM case. The general trend is for <!>_ to become larger on 
small scales. We also notice that the QS approximation seems to provide a very good match 
to the simulation results for A on all scales. Finally, as in the self-accelerating case, we see 
that the QS approximation provides reasonably accurate results (with errors < 5%) on scales 
k > OMhMpc-K 
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